Entropy Maximization and Instability of Uniformly Magnetized Plasma 
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The regime where a uniformly magnetized plasma is unstable to the spatial perturbation in the 
magnetic field is explored. In this regime, the entropy of the plasma system is not maximized in 
the uniformly magnetized state. The physical implication is discussed in the context of the current 
generation, the magnetic reconnection, and the dynamo effect. 
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It is often assumed that the uniformly magnetized plas- 
mas are dynamically stable, which is not vulnerable to 
a localized spatial perturbation in the magnetic field. If 
this is not the case, our premise on the uniform magnetic 
field that is often taken for granted should be reconsid- 
ered carefully. 

The second law of the thermodynamics states that 
the equilibrium state of a Hamiltonian system maxi- 
mizes the entropy for a given volume and the total en- 
ergy. The equilibrium configuration of a magnetized 
one-component electron plasma can be determined by 
this principle. More specifically, the electron distribu- 
tion function can be obtained by maximizing the Shan- 
non entropy in the framework of the variational princi- 
ple, under the constraints of the system-wide conserved 
quantities such as the gross number of the electrons, the 
volume, the energy, and the magnetic moment. Consider 
two kinds of plasma, one is with a uniform magnetic field 
Bi = Bqz and the other with a spatially-varying field of 
B2 = Bq(1 + ft cos(kx))z, where f3 is a small number. 
If the entropy of the latter system is larger than that of 
the former, it would attest that the plasma with a uni- 
form magnetic field could transition to a spatially varying 
state. The goal of this paper is to demonstrate that this 
is indeed the case in certain regime. 

Let us name the plasma with the uniform (spatially 
varying) magnetic field by the plasma U (V), and de- 
note the corresponding entropy by Su (Sy)- We assume 
that both plasmas have the same system-wide quantities, 
including the total energy, the magnetic moment, the to- 
tal number of the electrons and the volume. The total 
energy for the plasma U with a uniform magnetic field 
Bqz is 
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where m e is the electron mass and the summation of done 
over all the electrons. The interaction between the elec- 
tron magnetic moment and the magnetic field is ignored. 



The total magnetic moment is 

Mu = Ynr^ (2) 
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and the Shannon entropy is 

Su = -^2p q log(p q ), (3) 

k 

where c is the speed of light, p q is the probability that the 
plasma is in the state q, and the summation is done over 
all the possible states. Using the variation principle and 
the Lagrange multipliers l a and lb, 6 S + l a S Ejj + IbS Mjj = 
0, we obtain for each q 

6 Pq (- log(p g ) + l a E v + l b Mjj) = 0, (4) 

which leads to an anisotropic Maxwellian distribution 
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where vui_ {vu\\) is t ne perpendicular (parallel; z- 
direction) thermal velocity and no is the electron density. 
The temperature in each direction can be determined by 
the energy and the magnetic moment. 

Let us now consider the plasma V, with the magnetic 
field of £?o(l + (3 cos(kx))z where j3 <C 1. Following the 
same steps, the electron distribution is 
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where vy± C^Vll) is the perpendicular (parallel) thermal 
velocity 



v v± (xY 



V 



V± 



bp cos(kx) 
1+P cos(kx) 

2 



(7) 



°ux/ v 'u\\- Three unknowns, 
vv±i Vy\\ and ny, can be determined as a function of 



where 6 = 7 — ! with 7 = u^x/ u a|| 
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v u±, U(7||i and /3, given the constraints on the total 
number of electrons, energy and magnetic moment. The 
invariance in the number of electrons J fyd 3 v = J fud , v 
leads to 
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where (3 2 — ((3 2 cos(kx) 2 } = /3 2 /2, and the constraint on 
the total magnetic moment 



leads to 
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Lastly, the constraint on the energy is given as 
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J fv l m y = ^ + j fv l, 



This leads to the relationship 
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where A = (i?o/87r)/(m e w 2 / ||/2) is the ratio between the 
magnetic energy density and the electron parallel kinetic 
energy density. 

Note that, for a Maxwellian plasma, the Shannon en- 
tropy is proportional to 

S S - \og(n e /v 2 p v z ) = - log(n e )+log(t; 2 )+log(^), (13) 

where v 2 = v 2 + ia 2 . Defining Slog(n e ) = log(ny) — 
log(no), 6\og(v 2 ) = log(v vX ) - log(v% x ) and Slog(v z ) = 
log(wy||) — log^yy), we obtain from Eqs. dHJ), (fT0|) and 



51og(n e ) = 6(6 + 1)^ 2 , 

51og(i; 2 ) = -(26 2 + 36+l)/3 2 , 
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Finally, the entropy difference AS — Sv ~ Su is given 
by AS = (-(6 + l) 2 - A/2 + 7 (26 + 1)) /3 2 . From the 
relationship 7 = 6+1, this can be further simplified to 



AS = 7 2 -7- A/2, 
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which is the major result that our argument is based on. 

One important constraint that needs to be consid- 
ered is the self-consistency of the spatially varying mag- 
netic field BqP cos(kx) and the current generated by 



the magnetic moment j m = V x m, where m — 
(J fv(i r ncv 2 /B(x))d 3 v) z is the magnetic moment den- 
sity. In order for the spatially varying magnetic moment 
to be generated by the current of the spatially varying 
magnetic field, the relationship <5B = (47r/c)m, origi- 
nated from the Maxwell equation V x <5B = (47r/c)j m , 
needs to be imposed. Equating the relationship to the 
first order in /3, we obtain 
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The left-hand side of the above equation is exactly 7 2 /A, 
and Eq. (fTo]) is nothing but A = 7 2 . The condition 
AS > is a sufficient condition for the instability of the 
plasma U. Using Eqs. (fT"5|) and (|T6| . we obtain 7 > 2, 
which is the same condition with the classical Wcibel 
instability [1 Q. In order for our analysis to be valid, 
the wave vector k should be less than the inverse of 
the typical electron gyro-radius so that kvu±/uj ce < 1, 
where uj ce is the gyro-angular frequency While the 
unstable regime identified here is the same as for the 
Weibel instability, the physical origin is different. In the 
Weibel instability, the range of the unstable wave vector 
is < k < y // "f 2 /2 — l(u; po /c); however it is much wider 
in our theory, < k < r~ , where r g is the gyro-radius. 
It would be interesting to examine whether the instability 
identified here is much more explosive than the Weibel 
instability. 

The discussion on the constraint A = 7 2 is in order. 
It is imposed on the condition that the spatially varying 
part of the magnetic field is originated from the current 
generated by the electron magnetic moment. While this 
is a necessary condition for a self-sustaining system, it 
might be possible to drive the current from outside using 
various methodologies 0, 0| without perturbing the total 
energy or the magnetic moment of the plasma so that the 
condition can be lifted from the constraint. If A is very 
close to zero, then the instability can exist for a nearly 
isotropic Maxwellian plasma since the plasma is unstable 
even when 7=1. 

It should be noted that as the magnetic moment is 
not conserved in the presence of the collisions, the result 
presented here is valid only in the time scale faster than 
the electron collision rate. In a collisionless system, the 
magnetic moment distribution, as well as the total mag- 
netic moment itself, is conserved. As the plasmas U and 
V do not have the same distribution, the possibility of 
the dynamic connectivity between U and V should be 
addressed. One important point to be noted is that due 
to the spatial variation of the magnetic field, an elec- 
tron's magnetic moment is ambiguous up to the order 
of r g k. As a simplified example, consider B(x) = B + z 
if x > and B(x) = B-z if as < 0. An electron with 
the gyro-center at x = has two different magnetic mo- 
ments when it is moving at x > and x < 0. Although 
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the total magnetic moment is conserved (in the statis- 
tical sense), the final detailed distribution of the mag- 
netic moment is path-dependent upon how the magnetic 
field configuration changes from B z to B(l+/3cos(kx))z. 
Consider the two cases where the magnetic field changes 
from B z to B (l + /3 2 cos(2/cx))z and then finally to 
Bq(1 + f3cos(2x))z, and where it changes from Bqz di- 
rectly to Bq(1 + f3cos(kx))z. The distributions of these 
two cases cannot be the same. In other words, there is 
an infinite number of degrees of freedom for paths among 
which we might be able to choose a path that connects 
the plasma U and plasma V. While the above argument 
does not guarantee the dynamic connectivity between the 
plasma U and the plasma V, it shows that the chance 
that there exists a path from the plasma U to the plasma 
V, where the plasma V is very close to the plasma V, 
is high. 

With the caution on the dynamical connectivity be- 
tween the plasmas U and V, there could be implications 
of our theory on the dense plasma. As the plasma gets 
squeezed or expanded as in the inertial confinement fu- 
sion process [j| , the z-pinch plasma or solar corona mag- 
netic reconnection region Q, the plasma may cross the 
boundary given in Eq. ([T6|) with 7 > 2. If it does occur, 
as it crosses the boundary, the plasma would undergo 
one of the following three transitions before reaching to 
a stable state. First, the plasma may develop a spatially- 
varying magnetic field, as considered in this paper. This 
is possible only with the spatially- varying current genera- 
tion. Second, the plasma may convert the magnetic field 
energy into the electron kinetic energy, as in the case of 



the magnetic reconnection. Finally, the plasma may con- 
vert the electron kinetic energy into the magnetic field 
energy, as in the case of the dynamo effect. Our analysis 
does not answer which transition should occur. 

When the magnetic field is highly intense, as in the case 
of the astrophysical plasma, the quantum Landau level 
and other quantum effects become significant 
such a case, the ground state may not be at its minimum 
(maximum) energy (entropy) level when the magnetic 
field is uniform. This could leads to the interesting phe- 
nomena. Complication would be the quantum diffraction 
and degeneracy [1, 0-Q • 

We would like to thank Prof. Fisch on the various 
helpful discussion on the subject. 
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